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We illuminate the importance of a self-consistent many-body treatment in calculations of vac- 
uum polarization corrections to the energies of atomic orbitals in many-electron atoms. Including 
vacuum polarization in the atomic Hamiltonian causes a substantial re-adjustment (relaxation) of 
the electrostatic self-consistent field. The induced change in the electrostatic energies is substantial 
for states with the orbital angular momentum I > 0. For such orbitals, the relaxation mechanism 
determines the sign and even the order of magnitude of the total vacuum polarization correction. 
This relaxation mechanism is illustrated with numerical results for the Cs atom. 



PACS numbers: 31.30.Jv, 31.15.Ne, 31.15.Md, 31.25.Eb 



Compared to hydrogenic one-electron systems, calcu- 
lation of radiative corrections for many-electron atoms 
brings in an additional layer of complexity: a strong 
Coulomb repulsion between the electrons. The problem 
is especially challenging for neutral many-electron atoms, 
where the interaction of an outer-shell electron with other 
electrons is comparable to its interaction with the nu- 
cleus. At the same time, a reliable calculation of radia- 
tive corrections for a heavy neutral system is required 
in evaluation of the parity non-conserving (PNC) ampli- 
tude in the 55-electron 133 Cs atom. Here it has been 
only recently realized that the sizes of radiative correc- 
tions P, QJd) QL 01 are comparable to the experimental 
error bar jg of 0.35% and, together with the Breit cor- 
rection 0, dramatically affect agreement (or disagree- 
ment with the Standard Model of elementary parti- 
cles. 

A systematic approach to the problem of radiative cor- 
rections in strongly correlated systems is to start from a 
Furry representation based on a self-consistent electronic 
potential [9j- This potential takes into account the fact 
that an electron moves in an average field created by 
both the nucleus and other electrons. Based on this idea, 
a program of calculating radiative corrections to PNC 
amplitudes have been put forth by Sapirstein et al. p|. 
Kuchiev and Flambaum Q and Milstein et al. Q pur- 
sue a more qualitative approach using an independent- 
electron approximation. We believe that the question 
of an interplay between correlations and radiative cor- 
rections is yet to be addressed. While here we do not 
compute the PNC corrections, we illuminate a situation 
where disregarding correlations would lead to a substan- 
tial error in determining radiative correction: a radiative 
correction changes sign and even the order of magnitude 
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when the presence of other electrons is accounted for. 

In particular, we consider vacuum polarization (VP) 
corrections to energies of atomic states. To the leading 
order in aZ the VP may be accounted for by introducing 
the Uehling potential Uyp(r) into the atomic Hamilto- 
nian. This potential is attractive, and for a hydrogen- 
like ion the resulting VP corrections to the energies are 
always negative. For a complex atom, we find by con- 
trast that, for orbitals with I > 0, the total correction 
is positive. Briefly, the reason for such a counterintu- 
itive effect is due to a readjustment of atomic orbitals 
when the Uyp (r) potential is added to the self-consistent 
Dirac-Hartree-Fock (DHF) equations. The innermost Is 
orbitals are "pulled in" by the short-ranged VP poten- 
tial, leading to a decrease of the effective nuclear charge 
seen by the outer orbitals and thus to an increase of the 
electrostatic energy of these orbitals. Since for orbitals 
with I > 0, overlap with Uyp(r) and thus the lowest 
order correction are small, the resulting indirect "relax- 
ation" contribution dominates the total VP correction to 
the energies. In the following we will present numerical 
results supporting this relaxation mechanism. Atomic 
units (H = \e\ — m e = 1) are used throughout. 

Because of our interest in PNC in Cs, below we illus- 
trate the relaxation effect with numerical results for this 
atom; however, the relaxation mechanism is also appli- 
cable in the cases of other many-electron atoms. We also 
notice that the relaxation mechanism described here is 
similar to that observed in calculations of the Breit cor- 
rections [IbllTij. 

The conventional many-electron Hamiltonian may be 
represented as 



(i) 



where the single-particle Dirac Hamiltonian is 

h (i) = c(cti -Pi) + (3iC 2 + Kuc(^) • 

The nuclear potential V nuc (r) is obtained from the nu- 
clear charge distribution p nuc (r); which is we approxi- 
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mate by the Fermi distribution 

Pnuc( r ) = Z j 77 \ / -i ; (3) 

1 + exp[(r — c)/a\ 

where po is the normalization constant, c and a are the 
nuclear parameters. In our the numerical example for 
133 Cs, we use c = 5.6748 fm and a = 0.52 fm. 

A common starting point for describing a multi- 
electron atom is the self-consistent field method. Here 
the many-body wave-function is approximated by a 
Slater determinant constructed from single-particle or- 
bitals (bi-spinors) Uk(r). The orbitals are obtained by 
solving self-consistently the eigenvalue equations 

(h + E/dhf) Wfc(r) = e k u k (r) , (4) 

where £/dhf is the traditional DHF potential which de- 
pends on the orbitals occupied in the Slater determinant. 
The DHF energies for the core and several valence or- 
bitals of Cs are listed in Table [I] 

TABLE I: Vacuum polarization corrections to binding ener- 
gies in neutral Cs (Z = 55). Here e n i, are the DHF en- 
ergies, are the expectation values of the Uehling po- 
tential (Eq.©), and &°; HF are the VP corrections with the 
correlations included (Eq.JHJ). All quantities are given in 
atomic units, 1 a.u. = 27.21138 eV, and notation x[y] stands 
for x x 10 w . 



Orbital 




Se {1) 










core 


orbitals 








lsi/2 


-1330.396958 


-2.853 


-1] 


-2.782 


-1] 


2si/ 2 


-212.597116 


-3.392 


-2] 


-3.267 


-2] 


2pi/ 2 


-199.428898 


-1.510 


-3] 


5.406 


-4] 


2^3/2 


-186.434858 


-1.650 


-4] 


1.690 


-3] 


3«l/2 


-45.976320 


-6.868 


-3] 


-6.581 


-3] 


3pi/ 2 


-40.448097 


-3.339 


-4] 


1.987 


-4] 


3P3/2 


-37.893840 


-3.719 


-5] 


4.609 


-4] 


3d 3 / 2 


-28.309043 


-1.839 


-7] 


4.531 


-4] 


3d 5 / 2 


-27.774710 


-4.370 


-8] 


4.425 


-4] 


4si /2 


-9.514218 


-1.457 


-3] 


-1.397 


-3] 


4pi/2 


-7.446203 


-6.726 


-5] 


8.097 


-5] 


4P3/2 


-6.920865 


-7.506 


-6] 


1.355 


-4] 


4d 3/2 


-3.485503 


-3.440 


-8] 


1.153 


-4] 


4d 5/2 


-3.396788 


-8.100 


-9] 


1.129 


-4] 


5S1/2 


-1.490011 


-2.057 


-4] 


-2.050 


-4] 


5pi/ 2 


-0.907878 


-7.773 


-6] 


2.035 


-5] 


5^3/2 


-0.840312 


-8.395 


-7] 


2.757 


-5] 


valence states 








6S1/2 


-0.127380 


-1.054 


-5] 


-1.159 


-5] 


6^1/2 


-0.085616 


-1.942 


-7] 


2.284 


-7] 


6P3/2 


-0.083785 


-2.180 


-8] 


4.513 


-7] 


7S1/2 


-0.055190 


-2.896 


-6] 


-3.143 


-6] 


7P1/2 


-0.042021 


-6.957 


-8] 


8.150 


-8] 


7P3/2 


-0.041368 


-7.873 


-9] 


1.606 


-7] 



The polarization of the vacuum by the nucleus mod- 
ifies the nuclear electric field seen by the electrons. To 
the leading order in aZ, the VP may be conveniently de- 
scribed with the Uehling potential, which for a point-like 



nucleus of charge Z reads 

rrD c , , 2 aZ f°° , rr, / 1 1 \ T 2r ' 

(5)" 

This potential must be folded with the nuclear charge 
distribution, 



Uy P (r) = 




We approximated p n uc(^) with the Fermi distribution, 
Eq. 0. In the numerical evaluation of the extended- 
nucleus Uehling potential, we employed the routine from 
Ref. 0. The Uehling potential £/yp (r) generated by 
the Cs nucleus is shown in Fig. ^ Notice that the actual 
range of this potential is a few nuclear radii (instead of 
Compton wavelength A e ~ 384 fm) , because the potential 
for a point-like charge, Eq.©, diverges logarithmically as 
r — > 0; therefore the folded potential Uyp is dominated 
by the contributions accumulated inside the nucleus. 

o.oo 
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FIG. 1: Uehling potential for 133 Cs. Notice that the radius 
of the innermost Is orbital is about 10 3 fm, much larger than 
the effective range of the VP potential. 

How does one compute the VP corrections Set to the 
energies of the atomic orbitals? Below we consider two 
possibilities: (i) lowest-order perturbative treatment, 

fcW = (u k \Uvp\u k ), (6) 

and Qi) the self-consistent approach. Indeed, as in 
Ref. [3], the VP potential may be introduced into the 
DHF equations, 

(h + U VP + C^ HF ) u' k (r) = e' k u' k {r) , (7) 

and a set of new energies e' k and orbitals u' k (r) is ob- 
tained. Notice that the DHF potential is modified as 
well, since it depends on the new set of the occupied or- 
bitals u' k (r). The correlated VP correction to the energy 
of the orbital k is simply 

fe? HF = 4 - eu ■ (8) 
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Additionally, we carried out an independent corre- 
lated calculation in the framework of the linearized Cou- 
pled DHF approximation 0] , which is equivalent to the 
random-phase approximation (RPA). This approxima- 
tion describes a linear response of the atomic orbitals 
to the perturbing interaction, i.e. the VP potential. Nu- 
merical values obtained from the linearized coupled DHF 
calculations were in close agreement with the full DHF 
results. 

The numerical results of our calculations are presented 
in Table [I] While analyzing this Table, we observe that 
the lowest order corrections, Se^ , are always negative, 
reflecting the fact that the Uehling potential is attractive 
(see Fig. H). Owing to the short-ranged nature of VP, 
and the fact that only the s-orbitals have a significant 
overlap with the nucleus, the corrections to the energies 
of I = orbitals are much larger than those for I > 
orbitals. As to the correlated corrections, they differ 
quite substantially from the lowest order-corrections. A 
comparison of Eq. Q and Eq. Q reveals the origin of 
this discrepancy: the perturbation, in addition to the 
Uehling potential, contains a difference between the two 
DHF potentials 



5U = U V p + (U[ 



DHF 



Udhf) 



(9) 



For orbitals with I > 0, where the first term above is 
small, the modification of the DHF potential contributes 
significantly to the VP energy corrections. 

The modification of the DHF potential induced by the 
vacuum polarization is clearly a many-body effect, not 
present in hydrogen-like system. Such an effect has been 
explored before,for example in calculations of the Breit 
corrections [jjj, uM > an d h is commonly referred to as a 
relaxation mechanism. Let us illustrate this relaxation 
mechanism. Denoting the correction to the occupied or- 
bital wave functions as Xa(f) = u ' a ( r ) ~ u a(f), we write 

(£W - Uuhf) (r) « 

E / Xl(r')r-^u a (r')dr' + 
„ J \r-r'\ 

E / u a( r ') u ~ r n Xa(r')dr' - exchange , 

where we discarded contributions non- linear in Xa( r )> 
and "exchange" denotes non-local part of the perturba- 
tion. The first two (direct) terms can be interpreted as an 
electrostatic potential produced by a perturbation 6p c \(r) 
in the radial electronic density 



Pei(r) 



1 



47rr 2 ^ 



u\{r)u a (r) . 



We plot both the electronic density p e i(r) and the VP- 
induced perturbation 5p e \(r) in Fig. [5J The minima of 
p c \(r) correspond to positions of the electronic shells, 
marked on the plot by their values of principal quantum 
number n. 



The figure |3 may be interpreted in the following way: 
the s orbitals are "pulled in" by the attractive Uehling 
potential closer to the nucleus. As a result, screening of 
the nuclear charge by the inner orbitals becomes more 
efficient. For example, the modification of the effective 
charge felt by the n = 2 electrons is simply the area un- 
der the 5p e \(r) curve, accumulated between r = and 
the radius of the shell (r w 0.08ao); from Fig.|21it is clear 
that the induced modification of the effective charge for 
the n = 2 shell has a negative sign. Such an enhanced 
screening leads to a reduced attraction of the electrons by 
the nucleus and to the increase in the energy of the outer 
electrons. From Table H] we see that this indirect relax- 
ation contribution to the energy may be well comparable 
to the direct VP correction, 8e£ . While for I = orbitals 
the direct correction gives a reasonable estimate, for all 
orbitals with / > 0, the neglect of the relaxation would 
lead to even qualitatively incorrect result. Moreover, the 
higher the orbital angular momentum, the smaller is the 
direct correction, and the more important is the relax- 
ation mechanism. For example, for 4c? orbitals the VP 
correction in the lowest order is four orders of magnitude 
smaller than the correlated result. 
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FIG. 2: Perturbation of the electronic radial charge distribu- 
tion 6p e \(r) (solid line) for Cs atom due to vacuum polariza- 
tion by the nucleus. We also show the unperturbed density 
p e i(f) multiplied by a factor of 10 -3 (dashed line). The min- 
ima of p e i(r) correspond to positions of the electronic shells, 
marked on the plot by their values of the principal quantum 
number n. 



To summarize, here we illuminated the importance 
of the self-consistent many-body treatment in calcu- 
lations of vacuum polarization corrections. Including 
the VP Uehling potential into the atomic Hamiltonian 
causes re-adjustment (relaxation) of the electrostatic self- 
consistent field. The induced change in the electrostatic 
energies is substantial for states with the orbital angular 
momentum / > 0. As illustrated in our numerical results 
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for Cs, the relaxation mechanism determines the sign and work of W.R.J, was supported in part National Science 
even the order of magnitude of the total VP correction Foundation Grant No. PHY-0 1-39928. 
for orbit als with I > 0. 
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